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We give a review of modern theoretical understanding of the problem of thermal phase 

■ transition in QCD. The existence or non-existence of such transition depends on the nature 
CN ■ of the gauge group, the number of light quark flavors , and on the value of quark masses. 

■ Numerical lattice measurements indicate that the phase transition does not occur for the 
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experimentally observed values of quark masses. 



The properties of QCD medium at finite temperature have been the subject of intense study 
during the last 15 years. It was realized that the properties of the medium undergo a drastic 
change as the temperature increases. At low temperatures, the system presents a gas of 
colorless hadron states — the eigenstates of the QCD hamiltonian at zero temperature. When 



the temperature is small, this gas is composed mainly of pions — other mesons and baryons 
have higher mass and their admixture in the medium is exponentially small ~ exp{— M/T}. 
At small temperature, also the pion density is small — the gas is rarefied and pions practically 
do not interact with each other. 

However, when the temperature increases, pion density grows, the interaction becomes 
strong, and also other strongly interacting hadrons appear in the medium. For temperatures 
of order T ~ 150 Mev and higher, the interaction becomes so strong that the hadron states 
do not present a convenient basis to describes the properties of the medium anymore, and no 
analytic calculation is possible. 

On the other hand, when the temperature is very high, much higher than the characteristic 
hadron scale Hhadr ~ 0-5 Gev, theoretical analysis becomes possible again. Only in this range, 
the proper basis are not hadron states but quarks and gluons — the elementary fields entering 
the QCD lagrangian. For high temperatures, a characteristic energy of quarks and gluons 
travelling through the medium is also high, the effective coupling constant is small, and the 
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system presents the quark- gluon plasma. Its properties are in many respects very similar 
to the properties of the usual non-relativistic plasma involving charged particles with weak 
Coulomb interaction. The only difference is that quarks and gluons carry not the electric, 
but color charge. The perturbation theory in the coupling constant can be developed and 
many thermodynamic (such as free energy) and kinetic (such as viscosity) characteristics of 
the medium can be analytically evaluated 

Thus, the properties of the system at low and at high temperatures have nothing in com- 
mon. A natural question arises: What is the nature of the transition from low-temperature 
hadron gas to high- temperature quark- gluon plasma ? Is it a phase transition ? If yes, what 
is its order ? I want to emphasize that this question is highly non-trivial. A drastic change in 
the properties of the system in a certain temperature range does not guarantee the presence of 
the phase transition point where free energy of the system or its specific heat is discontinuous. 
Recall that there is no phase transition between ordinary gas and ordinary plasma. 

We shall see in the following that, as far as the real QCD with particular values of quark 
masses is concerned, the answer is probably negative. What really happens is not the phase 
transition but a sharp crossover — "almost" a second-order phase transition. However, the 
real phase transition does occur in some relative theories — in pure Yang- Mills theory (when 
the quark masses are sent to infinity) and in QCD with 2 or 3 exactly massless quark flavors. 

There are at least 4 reasons why this question is interesting to study: 

1. It is just an amusing theoretical question. 

2. Theoretic conclusions can be checked in lattice numerical experiments. Scores of papers 
devoted to lattice study of thermal properties of QCD have been published. 

3. Perhaps, a direct experimental study would be possible on RHIC — high-energy ion 
collider which is now under construction. I'll discuss the possibility to observe a beautiful 
effect, the so called disoriented chiral condensate in the end of the lecture. 

4. During the first second of its evolution, our Universe passed through the stage of high-T 
quark-gluon plasma which later cooled down to hadron gas (and eventually to dust and 
stars, of course). It is essential to understand whether the phase transition did occur 
at that time. A profound first-order phase transition would lead to observable effects. 
We know (or almost know — the discussion of this question has not yet completely died 
away) that there were no such transition. But it is important to understand why. 

Note that there is also a related but different question — what are the properties of relatively 
cold but very dense matter and whether there is a phase transition when the chemical potential 
corresponding to the baryon charge rather than the temperature is increased. This lecture will 
be devoted exclusively to the thermal properties of QCD, and we shall assume zero baryon 
charge density. 
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2 Pure Yang-Mills theory: deconfinement phase tran- 
sition. 



This is the system where the phase transition from the glueball phase to the gluon plasma 
phase does occur. This result has been obtained long ago by Polyakov 0] and Susskind Q. 
On the heuristic level, the reasoning is the following: 

We know (for real QCD — from experiment, and for pure YM theory — from theoretical 
arguments and from lattice measurements) that the theory enjoys confinement at low tem- 
perature. That means that the potential between the test heavy quark and antiquark grows 
linearly at large distances: 

T = 0: V QQ (r)~o-r, r ^ oo (2.1) 

On the other hand, at high temperature when the system presents a weakly interacting plasma 
of gluons, the behavior of the potential is quite different: 

T » Vhadr : V QQ (r) ~ l^l e ~m D r (2 . 2) 

Here ran ~ gT is the Debye mass, and the potential is the Debye screened potential much 
similar to the usual Debye potential between static quarks in non-relativistic plasma. There 
is no confinement at large T. There should be some point T c (the critical temperature) where 
the large r asymptotics of the potential changes and the phase transition from the confinement 
phase to the Debye screening phase occurs. 

These simple arguments can be formulated in a rigorous way. Consider the partition 
function of the system written as the Euclidean path integral. It is known since Matsubara 
that, at finite T, the fields are defined on the cylinder: Euclidean time r lies within the range 
< r < (3 = 1/T , and one should impose periodic boundary conditions on the gluon fields: 

A»(x,/J)=A£(x,0) (2.3) 

Let us choose a gauge where Aq is time- independent. Introduce the quantity called the 
Polyakov loop 

P(x) = i-Tr exp{ig(3A a (x)t°} (2.4) 
It is just a Wilson loop on the contour which winds around the cylinder. Consider the correlator 

C T (x)=<P(x)P*(0)> T (2.5) 



One can show 0] that the correlator (2.5) is related to the free energy of the test heavy 



quark-antiquark pair immersed in the plasma. 



Ct(x) = 1<*P{-0F® (r)} + iexp{-/JF$( r ) } (2.6) 
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where r = |x|. Fgl(r) and Fgl(r) are free energies of test quark-antiquark pairs (alias static 
potentials) in the triplet and, correspondingly, the singlet net color state. Let us take now the 
limit r — > oo. The quantity 

Ct(oo) = lim C T (r) (2.7) 

r — >oc 

plays the role of the order parameter of the deconfinement phase transition. At small T, FnQ(r) 
grow linearly at r — > oo and Ct(oo) = 0. At large T, free energies do not grow and Ct(oo) 
is some non-zero constant (if one would naively substitute in Eq.(|2.6|) the Debye form of the 



potentials (22), one would get CT> /Jhadr (oo) = 1, but it is not quite true because FnQ(r) 
involve also a constant depending on the ultraviolet cutoff of the theory. See H ||] for detailed 
discussion). There is a phase transition in between. 

What are the properties of this phase transition ? There are not quite rigorous but sugges- 
tive theoretical arguments based on the notion of "universality class" [|?J which predict different 
properties for different gauge groups. The main idea is that the pure YM theory based on 
SU(2) color group has some common features with the Ising model (with global symmetry 
Z2), the theory with SU(3) gauge group — with a generalized Ising model (the Potts model) 
with the global symmetry Z3 etc. The Ising model has the second order phase transition, and 
the same should be true for pure SU(2) gauge theory. Systems with Zn symmetry display, 
however, the first order phase transition, and the same should be true for pure SU(N > 3) 
theory. The lattice data Q are in a nice agreement with this prediction. Also, critical indices 
of the second order phase transition were measured. Their numerical values are close to the 
numerical values of critical indices in the Ising model. 

2.1 Bubble confusion. 

There was a long-standing confusion concerning the nature of deconfinement phase transition 
in pure YM theory. It has been clarified only recently and I want to dwell on this question in 
more details. 

In scores of papers published since 1978, it was explicitly or implicitly assumed that one 



can use the cluster decomposition for the correlator (2.5) at large T and attribute the meaning 
to the temperature average < P >t- Under this assumption, the phase of this average can 
acquire N c different values: < P >t= C exp{27rik / N c } , k = 0, 1, . . . , N c — 1 which would 
correspond to N c distinct physical phases and to the spontaneous breaking of the discrete 
^TV-symmetry. In recent ||, the surface energy density of the domain walls separating these 
phases has been evaluated. 

However, the standard interpretation is wrong. In particular: 

1. Only the correlator (|2.5| ) has the physical meaning. The phase of the expectation value 
< P >t is not a physically measurable quantity. There is only one physical phase in the 
hot YM system. 
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2. The "walls" found in || should not be interpreted as physical objects living in Minkowski 
space, but rather as Euclidean field configurations, kind of "planar instantons" appearing 
due to non-trivial tti(G) = Zj^ where Q = SU(N)/Zn is the true gauge symmetry group 
of the pure Yang-Mills system. 

3. The whole bunch of arguments which is usually applied to non-abelian gauge theories 
can be transferred with a little change to hot QED. The latter also involves planar 
instantons appearing due to non- trivial 7ri[t7(l)] = Z. These instantons should not be 
interpreted as Minkowski space walls. 

It is impossible to present an adequate discussion of this issue in this short lecture. We 
refer the reader to || where such discussion was given. Here we restrict ourselves by outlining 
some heuristic physical arguments. 

Right from the beginning, one meets a puzzle. In standard approach, Zn is broken spon- 
taneously at high temperature and restored at low temperature. This is very strange and 
unusual. The opposite is much more common in physics. There are some models where sym- 
metry breaking survives and can even be induced at high temperature [[h]], but the mechanism 
of this breaking is quite different from what can possibly occur in the pure YM theory. 

The second observation is that Zn symmetry which is presumably broken in the decon- 
finement phase is just not there in the continuum theory: the gluon fields are not transformed 
under the action of the center. Z^ - symmetry is present in the standard lattice version of the 
theory but is absent, again, in the lattice theory involving adjoint matrices O ab = Tr{t a Ut b U^} 
rather than the unitary matrices U. 

It was mentioned earlier that the Polyakov loop expectation value < P >t as such has 
no physical meaning. Let us explain why. Being taken at face value, < P >t would measure 
the free energy of a single fundamental heavy source immersed in the system |ll]]: < P >t= 
exp{— PFt}- As a matter of fact, nonzero phase of < P >t would correspond to the complex 
free energy which is an obvious nonsense. But the point is that one just cannot put a single 



fundamental source into the system due to the Gauss law constraint [12]. The net color charge 
should be zero, and a fundamental source cannot be screened by gluons — the only dynamic 
fields in the lagrangian of the theory and in the heat bath f^] 

What one can well do is to immerse a heavy quark-antiquark pair and measure thereby 
the correlator ( |2.5D which is physical. Or, say, for the SU(3) gauge group, one can immerse 3 
heavy quarks at different points and measure < P(x)P(y)P(0) >t which is physical, again, 
but, in contrast to < P >t, does not involve the phase uncertainty. 

Actually, the delusion of spontaneous Zn breaking in pure YM theory persisted for so long 



3 To be quite precise, the net color charge can be made non-zero due to a boundary term at spatial infinity. In a 
finite spatial box, such a term can appear when non-standard (non-periodic) boundary conditions are chosen. But 
then the phase of < P >t would be exactly determined by these boundary conditions. As the physical properties of 
the theory cannot depend on b.c. when the box is large enough, it is just another way to say that the phase is not 
physical. 
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because people habitually described this system in terms of ^4q( x ) ( more exactly — in terms 
of f2 x = exp{i/3gAg(x)i a }) which are not the dynamic variables entering the hamiltonian but 
the variables dual to the Gauss law constraints. A close analogy can be drawn with the Ising 
model in 2 dimensions. When expressed in terms of the original spin variables cjj, the system is 
ordered at low temperatures and disordered at high temperatures — the spontaneously broken 
Z2 - symmetry is restored there. One can make, however, the Kramers- Wannier transforma- 
tion and describe the system in terms of the dual or "disorder" variables rji When the 
normal temperature is high, the dual temperature is low, and the dual hamiltonian H*[r]i] 
describes, indeed, the system where Z2 symmetry is broken at high temperatures and restored 
at low temperatures. But the variables rji are not physical observables and, as far as any 
Gedanken physical experiment is concerned, Z2 - symmetry in the Ising model is restored at 
high temperatures. 

To summarize, there is only one physical phase at high T. Its properties are relatively 
simple — it is the weakly interacting plasma of gluons. The description in terms of dual 
variables is useful for some purposes (e.g. the universality class arguments of Ref. [0 which 
predict the order of the phase transition are based on the dual description), but one should be 
very careful not to read out in it something which is not in Nature. 

3 QCD with massless quarks. 

If the theory involves besides gluons also quarks with finite mass, the static interquark potential 
Vqq{t) does not grow at large distances anymore even at T = 0. Dynamic quarks screen the 
potential of static sources. One can visualize this screening thinking of the color gluon tube 
stretched between two static fundamental sources being torn apart in the middle with the 
formation of an extra quark-antiquark pair. Thus, in QCD with quarks, the Wilson loop 
average has the perimeter rather than the area law The correlator of two Polyakov loops 
( |2.5| ) tends to a constant at large distances universally at low and at high temperature, and 
this correlator cannot play the role of the order parameter of phase transition. 

Still, the phase transition can occur and does occur in some versions of the theory. It is 
associated, however, not with change in behavior of the correlator fl2.5|), but with restoration 
of chiral symmetry which is spontaneously broken at zero temperature. 

Consider YM theory with SU(3) color group and involving Nf massless Dirac fermions in 
the fundamental representation of the group. The fermion part of the lagrangian is 

/ 

where T>^ = — igA®t a is the covariant derivative. The lagrangian ( |3.1| ) is invariant under 

4 That does not mean that there is no confinement — as earlier, only the colorless states are present in the physical 
spectrum. But the behavior of the Wilson loop is not a good signature of confinement anymore. 
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chiral transformations of fermion fields: 



Ql,r A LtR q L:R (3.2) 

where qL,R = 4(1 ± 7 5 )<7 is the flavor vector with Nf components and Al^r are two different 
U(Nf) matrices. Thus, the symmetry of the classical lagrangian is Ui,(Nf) ® Un(Nf). Not 
all Nother currents corresponding to this symmetry are conserved in the full quantum theory. 
It is well known that the divergence of the singlet axial current = J2f Qfl^l^Qf is non-zero 
due to anomaly: 

~ g 2 e^G%G a ap (3.3) 

Thus, the symmetry of quantum theory is SUL(Nf)®SUn(Nf)®Uv(l)- It is the experimental 
fact that (for Nf = 2,3, at least) this symmetry is broken spontaneously down to Uy{Nf). 
The order parameter of this breaking is the chiral quark condensate < J2fQfQf >o- This 
spontaneous breaking leads to appearance of the octet of pseudoscalar Goldstone states in the 
spectrum. Of course, in the real QCD the quarks are not exactly massless, the mass term 
is not invariant with respect to the symmetry (|3.2| ) but only under Uv(Nf). As a result, in 
real World we have the octet of light (but not massless) pseudo- Goldstone pseudoscalar states 
(tt, K, rj). But the small mass of pseudogoldstones and the large splitting between the massive 
states of opposite parity (p/A\, etc.) indicate beyond reasonable doubts that the exact chiral 
symmetry (|3.2|) would be broken spontaneously in the massless case. As the masses of the 



strange and, especially, of u- and d- quarks are small [14| , the mass term in the lagrangian 



can be treated as perturbation. E.g. the pion mass satisfies the relation 

Flm\ = (m u + m d )\ < uu > | (3.4) 

(.Ftt = 93 Mev is the pion decay constant) and turns to zero in the chiral limit m U) d — ► 0. 
It is noteworthy that the symmetry breaking pattern 

SU L (N f ) ® SUfi(Nf) SU v (N f ) (3.5) 

depends crucially on the assumption that the gauge group involves at least 3 colors. For SU (2) 
color group where quarks and antiquarks belong to the same representation (the fundamental 
representation of the SU(2) group is pseudoreal: 2 = 2), the symmetry group of the lagrangian 
Q3.1| ) is much higher. It is U(2Nj) and involves also mixing between quarks and antiquarks. 
Ua(X) ' P ar t of this symmetry is anomalous and the formation of chiral condensate breaks 
spontaneously the remaining SU{2Nf) down to a simplectic group |jTq] : 

SU(2N f ) -» Sp{2N f ) (3.6) 

As a result, 2A r | — Nf — 1 Goldstone bosons living on the coset space appear. For Nf = 2, 
we have not 3 as usual, but 5 "pions". This fact is important to understand for people who 
would wish to study numerically on lattices the spontaneous chiral symmetry breaking with 
SU(2) gauge group. 
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As far as the thermal properties of the theory are concerned, the point is that a spon- 
taneously broken symmetry must be restored under a sufficient heating. There should be a 
critical temperature T c above which the fermion condensate < qq >t is zero. This is the 
temperature of phase transition and < qq >t is the order parameter associated with the 
transition. 

Note that the phenomenon of spontaneous chiral symmetry breaking is specific for theories 
with several light quark flavors. In the theory with Nf = 1 , the non- anomalous part of the 
symmetry of the lagrangian is just £7y(l). It stays intact after adding the mass term and after 
taking into account the formation of the condensate < qq >. The condensate is still formed, 
but it does not correspond to spontaneous breaking of any symmetry and need not vanish at 
high temperature. So, it does not. At high temperatures when the effective coupling is small, 
it can be evaluated semiclassically in the instanton approach |l6|, and one can show that 
it falls down as a power of temperature and never reaches zero. 

There is no phase transition in QCD with only one light or massless flavor. But it does 
occur when the number of massless flavors Nf is 2 or more. 

The melting down of quark condensate can be studied analytically at low temperature 
when the medium presents a rarefied weakly interacting gas of pions with low energies. Their 
properties are described by the effective chiral lagrangian 

C = ^Tr{d,Ud,U^ + ... (3.7) 

where U is the SU(Nf) matrix and the dots stand for higher derivative terms and the terms 
involving quark masses. When the characteristic energy and the quark masses are small, the 
effects due to these terms are suppressed and a perturbation theory (the chiral perturbation 
theory |l8|]) can be developed. In JT^], the temperature dependence of < qq >t has been 
determined on the 3-loop level. In the approximation where only the presence of pions in the 
heat bath is taken into account and the effects due to non-zero m u and are neglected, the 
result has a rather simple form 



< qq >t=< qq >o 



rj>2 rp4 rp& 

■ In \~ ... 



8F2 384Fi T 



(3.8) 



The constant A depends on the higher-derivative terms in the effective lagrangian and can be 
fixed from experiment: A ~ 500 ± 100 Mev. The dependence fl3.8j ) together with the curves 
where only only the 1 loop correction oc T 2 and 2 loop correction oc T 4 are taken into account 
(please, do not put attention to the "technical" curve marked o° = 0) is drawn in Fig. 1 taken 
from Ref. [FlI 
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The expansion in the parameter ~ T 2 /8F 2 makes sense when this parameter is small, 
i.e. when T < 100 — 150 Mev. Strictly speaking, one cannot extrapolate the dependence 
(|3.8| ) for larger temperatures, especially having in mind that, at T > 150 Mev, the heat bath 
includes a considerables fraction of other than pion hadron states. But as we know anyhow 
that the phase transition with restoration of chiral symmetry should occur, the estimate of the 
phase transition temperature (i.e. the temperature when < qq >t hits zero) based on such an 
extrapolation is not altogether stupid. This estimate is 

T c ps 190 Mev (3.9) 

(A more accurate treatment which takes into account non-zero m Uj a and also the presence of 
other mesons in the heat bath gives practically the same estimate as these two effects push T c 
in opposite directions and practically cancel each other.) 

4 Properties of phase transition. The real World. 

Making the estimate ( |3.9| ), we tacitly assumed that the phase transition is of the second order: 
only the derivative d < qq >t /dT but not < qq >t is discontinuous at the phase transition 
point. Let us discuss the question whether this assumption is valid and under what conditions. 

On the theoretical side, the situation is similar to that in pure YM case: not rigorous 
but suggestive arguments exist indicating that the phase transition is of the second order for 
2 massless flavors. When Nf > 3, the phase transition is probably of the first order. The 
arguments are the following |2Q j: 

The starting point is the observation that, in theories involving scalar fields, phase tran- 
sition of the first order often occurs when the potential involves a cubic in fields term. One 
can recall in the first place a cubic Van-der-Vaals curve P(p, T) which describe the first order 
water <-> vapor phase transition. The simplest field theory example is the theory of real scalar 
field with the potential 

V(4>) = \(tf-v 2 ) 2 -fi ( f> 3 (4.1) 

Assume for simplicity [i <C \v. At T = 0, the potential has one global minimum at <f) ~ 
■u + 3^/8A. At non-zero temperature, the term ~ XT 2 (p 2 is added to the effective potential. At 
high temperature T>ti, the minimum occurs at <j> = 0. One can be easily convinced that a 
local minimum at <j) = appears at some temperature T* when the local minimum at positive 
<j) still exists. The latter disappears at some larger temperature T**. In a certain temperature 
range, two minima of the potential, the old and the new one, coexist, one being a metastable 
state with respect to the other. This is exactly the physical situation of the first order phase 
transition. 

Let us go back to QCD. A direct application of this reasoning is not possible because the 
QCD lagrangian does not involve scalar fields. The effective chiral lagrangian ( [3.7D is also of 
no immediate use because higher-derivative terms which stand for dots cannot be neglected 
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in the region close to critical temperature. Suppose, however, that in the region T ~ T c some 
other effective lagrangian in Ginzburg-Landau spirit can be written which depends on the 
composite colorless fields 

= QRfQLf (4.2) 
A general form of the effective potential which is invariant under SUL(Nf) (g) SUn(Nf) is 

V[$] ~ ^iTr-t^} + g 2 (Tr{^}) 2 + £ 3 Tr{$$ t $$ t } + # 4 (det $ + det $ f ) + . . . (4.3) 

(the coefficients may be smooth functions of T). Now look at the determinant term. For 
Nf = 2, it is quadratic in fields while, for Nf = 3, it is cubic in fields and the effective 



potential acquires the structure similar to Eq.(4.1) which is characteristic for the systems with 



first order phase transition. A more refined analysis [2(J shows that the first order phase 
transition is allowed also for Nf > 4, but not for Nf = 2 where the phase transition is of the 
second order. 

It is even possible to argue that, for Nf > 3, one has not one but two phase transitions. The 
argument is based on the exact relation for the spectral density of Euclidean Dirac operator 
p(\) in zero-temperature QCD with Nf massless flavours at small but nonzero A. It is possible 



to show that it involves a non- analytic term in A at Nf > 3 [21 



where E = | < qq >o |. At nonzero temperature, both p(0) (i.e. the chiral condensate) and 
the coefficient of |A| in p(X) are changed. Eventually (at T 3> Phadr) they should vanish. One 
can imagine a situation when the chiral condensate turns to zero before the slope does. Or 
other way round. 

If lattice people would eventually observe two temperature phase transitions in QCD with 3 
massless quark flavors, I would be happy, of course. However, a canonical viewpoint that there 
is only one transition but of the first order may also be true. Existing lattice measurements 
favour this possibility. 

Up to now, we discussed only pure YM theory and QCD with massless quarks. But the 



quarks have non-zero masses: m u ~ 4 Mev, ~ 7 Mev, and m s ~ 150 Mev [14]. The question 
arises whether the non-zero masses affect the conclusion on the existence or non-existence and 
the properties of the phase transition. 

The experimental (i.e. lattice) answer to this question appears to be positive |22]| . In Fig. 
2, a phase diagram of QCD with different values of quark masses m s and m u = is plotted. 
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Let us discuss different regions on this piot. When the quark masses are large, quarks 
effectively decouple and we have pure YM theory with SU(3) gauge group where the phase 
transition is of the first order. When all the quark masses are zero, the phase transition is also 
of the first order. When masses are shifted from zero a little bit, we still have a first order 
phase transition because a finite discontinuity in energy and other thermodynamic quantities 
cannot disappear at once when external parameters (the quark masses) are smoothly changed. 

But when all the masses are non-zero and neither are too small nor too large, phase 
transition is absent. Notice the bold vertical line on the left. When m u = = and m s is 
not too small, we have effectively the theory with two massless quarks and the phase transition 
if of the second order. The experimental values of quark masses (the dashed circle in Fig. 2) lie 
close to this line of second order phase transitions but in the region where no phase transition 
occurs. It is the experimental fact as measured in Ref. |22| j. 

This statement conforms nicely with a semi-phenomenological theoretical argument of ref. 



1 23 ] which displays that even if the first order phase transition occurs in QCD, it is rather 
weak. The argument is based on a generalized Clausius-Clapeyron relation. In college physics, 
it is the relation connecting the discontinuity in free energy at the first-order phase transition 
point with the sensitivity of the critical temperature to pressure. The Clausius-Clapeyron 
relation in QCD reads 

1 dT c 

disc < qq >t= t^t disc e (4-5) 

T c dm q 

where disc e is the latent heat. The derivative J^- can the estimated from theoretical and 
experimental information of how other essential properties of QCD depend on m q and from 
the calculation of T - dependence of condensate at low temperature in the framework of chiral 
perturbation theory (see Fig.l and the discussion thereof). The dependence on quark masses 
is not too weak. From that, assuming that the discontinuity in quark condensate is as large 
as < qq >o (which is not true, of course), we get an estimate 

disc e < 0.4 GeV/fm 3 

which is rather small compared to the characteristic free energy density of the system in the 
vicinity of T c ~ 190MeV. 

Thus, latent heat of the first order phase transition (assuming it is there) must be small 
which means that the phase transition is likely to disappear under a relatively small pertur- 
bation due to nonzero m s . 

The question is not yet completely resolved, and independent lattice measurements are 
highly desirable. Most probable is, however, that, when temperature is changed, hadron gas 
goes over to quark-gluon plasma and other way round without any phase transition. There is, 
however, a sharp crossover in a narrow temperature range which is similar in properties to a 
second-order phase transition (the "phase crossover" if you will) . 
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5 Instantons and percolation. 



In the analysis in previous two sections, we relied on the fact that chiral symmetry is broken 
at zero temperature. It is an experimental fact in real QCD, but it is important to understand 
from pure theoretical premises why it is broken and what is the mechanism of its restoration 
at higher temperatures. 

A completely satisfactory answer to this question has not yet been obtained. The problem 
is that QCD at zero temperature is a theory with strong coupling and it is very difficult 
(may be impossible) to study the structure of QCD vacuum state analytically. However, a 
rather appealing qualitative physical picture exists which is based on the model of instanton- 
antiinstanton liquid and on the analogy with the so called percolation phase transition in doped 



superconductors [24]. We refer the reader to the Shuryak's book for the detailed discussion 
and elucidate here only crucial points of the reasoning. 

The starting point is the famous Banks and Casher relation |25[] connecting quark conden- 
sate to the mean spectral density of Euclidean Dirac operator p(X) at A ~ 0. Let us explain 
how it is derived. Consider the Euclidean fermion Green's function < q(x)q(y) > in a par- 
ticular gauge field background. Introduce a finite Euclidean volume V to regularize theory in 
the infrared. Then the spectrum of massless Dirac operator is discrete and enjoys the chiral 
symmetry: for any eigenfunction ip n {x) satisfying the equation "ptpn = X n i^ n , the function 
ip n = ^ipn is also an eigenfunction with the eigenvalue A n = — X n . 

The idea is to use the spectral decomposition of the fermion Green's function with a small 
but non-zero quark mass 

<q(x)q(y)> = ^MiMM (5.1) 
^ iX n -m 

Set x = y and integrate over d A x. We have 

V < qq >= -m £ -^-3 (5.2) 

where the chiral symmetry of the spectrum has been used and the contribution of the zero 
modes A ra = has been neglected (it is justified when the volume V is large enough |2q]). 
Perform the averaging over gauge fields and take first the limit V — > 00 and then the limit 
m — > 0. The sum can be traded for the integral: 

< qq >= -m f -^—dX = --p(0) (5.3) 
J X z + m z tt 

The rightmost-hand-side of Eq,(|5.3|) is only the non-perturbative m- independent part of the 
condensate . There is also a perturbative ultraviolet-divergent piece oc mA^ iir which is pro- 
portional to the quark mass, is related to large eigenvalues A and is of no concern for us 
here. 

Thus, the non-perturbative part of the quark condensate which is the order parameter of 
the symmetry breaking is related to small eigenvalues of Euclidean Dirac operator. There 
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should be a lot of them — a characteristic spacing between levels is SX ~ 1/(| < qq > \V) 
which is much less than the characteristic spacing 5X ~ 1/L for free fermions. 

The question is what is the physical reason for these small eigenvalues to appear. As 
far as we know, the first pioneer paper where a mechanism for generating small eigenvalues 
was proposed is Ref. [^7]] where small eigenvalues appeared as zero modes of monopole-like 
gauge field configurations. The disadvantage of this model is that the monopole configurations 
are static whereas it is natural to expect that characteristic gauge fields contributing to the 
Euclidean path integral at T = are more or less symmetric in all four directions with no 
particular axis being singled out. The model of instanton-antiinstanton liquid formulated 
in p8| (see in particular Ref. [28b] where the mechanism for spontaneous chiral symmetry 



breaking was suggested) and developed later in [29] is much better in this respect. 

The basic assumption of the model is that a characteristic gauge field contributing in QCD 
path integral is a medium of instantons and instantons as shown in Fig. 3. It is not a "gas" of 
Callan, Dashen, and Gross [30] because the interaction between instantons and antiinstantons 
bringing about a short-range correlations between instanton positions and orientations cannot 
be neglected. A "liquid" is a more proper term. 




Fig. 3. Instanton-antiinstanton liquid. 



Each individual instanton and antiinstanton involves a fermion zero mode [31]. Assuming 
the constant density of quasi-particles oc n\ a d r , the total number of zero modes in the Euclidean 
volume V is N ~ V^tadr- However, these are not exact zero modes. They are shifted from zero 
due to interaction between instantons and antiinstantons (a nonzero overlap between individual 
instanton and antiinstanton zero modes. Assuming their uniform spreading in the range of 
eigenvalues A A ~ l^hadn 

the volume density of quasi-zero modes is p(0) ~ N/(VAX) ~ /J% adr . 
Due to Eq. (|5.3[) , a non-zero quark condensate appears 0. 



5 The assumption of quasi-uniform spreading of eigenvalues is not so innocent. It probably holds only in the 
theory with several light dynamical quarks, but not in the quenched theory (Nf — 0) where it is natural to expect a 
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This picture is rather similar to what happens in a doped superconductor with high enough 
doping. When a characteristic distance between individual atoms of the admixture is not large, 
the wave functions of outer electrons of these atoms overlap, and the electrons can jump from 
site to site. If the set of atoms of admixture with a noticeable overlap of wave functions 
forms a connected network in the space, the electrons can travel through this network at large 
distances and the specimen is a conductor. Note that it is not a standard metal mechanism of 
conductivity when the medium is a crystal, has the long-range order, and the electron wave 
functions are periodic Bloch waves. Here the distribution of the dope whose electrons are 
responsible for conductivity is stochastic and wave functions are complicated. The essential is 
that they are delocalized. 

Thus, one can say that the vacuum of QCD is the "conductor" in a certain sense. For 
sure, there is no conductivity of anything in usual Minkowski space-time. Only the Euclidean 
vacuum functional has "conducting" properties. In principle, one can introduce formally the 
fifth time and write an analog of Kubo formula for conductivity in QCD, but the physical 
meaning of this "conductivity" is not clear. It is sufficient to say that, in a characteristic 
Euclidean gauge field background, the eigenfunctions of Dirac operator corresponding to small 
eigenvalues are delocalized. 

What happens if we heat the system ? The effective coupling constant g 2 (T) decreases, 
the action of individual instantons S = 8tt 2 /g 2 (T) increases, and the density of quasi-particles 
oc exp{— S} decreases. 

Let us look first at the doped superconductor when we decrease the density of admixture. 
Below some critical density, the set of atoms with essential overlapping of wave functions does 
not form a connected network in 3-dimensional space anymore. The electrons can no longer 
travel far through this network, wave functions become localized, and the specimen is an 



insulator. This is called the percolation phase transition (see e.g. [33] for detailed discussion). 

Likewise, there is a critical temperature in QCD above which instantons and antiinstantons 
do not form anymore a connected cluster with an essential overlap of individual fermion zero 
modes [what overlap is "essential" and what is not is a numerical question. For condensed 
matter systems (but not for QCD in this context) a computer estimates for the critical ad- 
mixture density has been performed]. At high temperatures, few remaining quasi-particles 
tend to form "instanton-antiinstanton molecules" (See Fig. 4). The individual zero modes are 
not spread out uniformly in the range AA ~ fJ-hadr as is the case at zero temperature where 
instantons and antiinstantons form an infinite cluster, but are just shifted by the value ~ Hhadr 
due to interaction in individual molecules. Small eigenvalues in the spectrum of Dirac operator 



are absent and the fermion condensate is zero 24]. 



singular behaviour of the spectral density near zero: p(A) ~ 1/A so that the "fermion condensate" (i.e. the vacuum 
expectation value < qq >q where quark fields are treated as external sources) is infinite f32[|. 
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Fig. 4. Gas of instanton-antiinstanton molecules (high T). 

Of course, this picture is too heuristic and qualitative. A serious quantitative study of 
fermion eigenvalues and eigenfunctions at non-zero temperature, and especially in the region 
T ~ T c has not yet been done. It is the task (a very interesting and important one) for future 
explorers. 

It is worthwhile to emphasize once more that this scenario of percolation phase transition 
leading to the molecular high-temperature phase is expected to hold only at Nf > 2. For Nf = 
1 with arbitrary small but nonzero fermion mass, molecules get ionized and the "medium" 
presents a very dilute instanton-antiinstanton gas — the instanton density involves a product 
of two small factors: exp{— 87r 2 /<7 2 (T)} and the fermion mass m. Differentiating logZ over 
m and sending m to zero, one gets a small but non-zero quark condensate (l^, |17]]. Cf. the 
analogous situation in the Schwinger model 

6 Disoriented chiral condensate. 

When we talked in previous sections about "experimental" tests of theoretical predictions, 
we meant numerical lattice experiment. It is the unfortunate reality of our time that the 
feedback between theory and real laboratory experiment has drastically deteriorated: what is 
interesting from theoretical viewpoint cannot very often be measured in laboratory and what 
can be measured is not interesting. 

However, speaking of the particular problem of the phase transition in QCD associated 
with chiral symmetry restoration, an intriguing possibility exists that a direct experimental 
evidence for such a transition can be obtained at the high-energy heavy ion collider RHIC 
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which is now under construction. After a head-on collision of two energetic heavy nuclei, a 
high temperature hadron "soup" is created. 

We do not call this soup the quark-gluon plasma because, even at RHIC energies, the 
temperature would not be high enough to provide a sufficient smallness of the effective coupling 
g 2 (T) and to make the perturbation theory over this parameter meaningful. What is important, 
however, is that, at RHIC energies, the temperature of the soup would be well above the 
estimate ( |3.9|) for the phase transition temperature. The high-temperature state created in 
heavy nuclei collision would exist for a very short time after which it expands, is cooled down 
and decays eventually into mesons. 

Let us look in more details at the cooling stage. At high temperature, the fermion conden- 
sate is zero. Below phase transition, it is formed and breaks spontaneously chiral symmetry. 
This breaking means that the vacuum state is not invariant under the chiral transformations 
( |3.2| ) and a direction in isotopic space is distinguished. What particular direction — is a matter 
of chance. This direction is specified by the condensate matrix 

=< QLfQRf > (6-1) 
For simplicity, we have assumed up to now that 

= ( 6 - 2 ) 

, but any unitary matrix can be substituted for Sfy (of course, it can be brought back in the 
form 5ff by a chiral transformation ). In different regions of space, cooling occurs indepen- 
dently and directions of condensate are not correlated. As a result, domains with different 
directions of condensate shown in Fig. 5 are formed (cf. cooling down of a ferromagnetic below 
the Curie point). 

In our World, we do not observe any domains, however. The direction of the condensate 
in all spatial points is identical. This is a consequence of the fact that u- and d- quarks have 
non-zero masses which break chiral symmetry explicitly, the vacuum energy involves a term 

E vac ~ Ti{M^} + cx . ( 6 .3) 

and the only true vacuum state is (|6.2| ) (in the basis where the quark mass matrix Ai is 
diagonal). 

However, the masses of u- and d- quarks are rather small and one can expect that the 
domains with "wrong" direction of the condensate are sufficiently developed during the cooling 
stage before they eventually decay into true vacuum (|T^) with emission of pions. [] 

This is a crucial assumption. A theoretic estimate of the characteristic size of domains they 
reach before decaying is very difficult and there is no unique opinion on this question in the 

6 Fig.5 implies the existence of several domains and describes better the physical situation immediately after the 
"phase crossover" in early Universe. Probably, the size of the hot fireball produced in collision of two nuclei is too 
small and the cooling occurs too fast for several domains to be developed. The popular "baked Alaska" scenario |35| 
implies the formation of only one domain with (generally) wrong flavour orientation. 
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literature. But if this assumption is true, we can expect to observe a very beautiful effect pi]] . 
From the true vacuum viewpoint, a domain with disoriented ^ff is a classical object — kind of 
a "soliton" (parentheses are put because it is not stable) presenting a coherent superposition of 
many pions. The mass of this quasi-soliton is much larger than the pion mass. The existence 
of such multipion coherent states was discussed long ago in pioneer papers |3(| but not in 
relation with thermal phase transition. 

Eventually, these objects decay into pions. Some of the latter are neutral and some are 
charged. As all isotopic orientations of the condensate in the domains are equally probable, 
the average fractions of ir°, tt + , and 7r~ are equal: < f T o >=< f w ± >= | as is also the case 
for incoherent production of pions in, say, pp collisions where no thermalized high-T hadron 
soup is created. 

But the distribution P{f) over the fraction of, say, neutral pions is quite different in the case 
of incoherent and coherent production. In incoherent case, P(f) is a very narrow Poissonic 
distribution with the central value < f w o >= 1/3. The events with f n o = or with f n o = 1 
are highly unprobable: P(0) ~ -P(l) ~ exp{— CiV} where N 3> 1 is the total number of pions 
produced. 

For coherent production, the picture is quite different. £jy is proportional to a SU(2) 
matrix. Factorizing over £7(1), one can define a unit vector in isotopic space E S 2 . The fraction 
of 7T° produced would be just / = cos 2 9 where 6 is a polar angle on S 2 . The probability to 
have a particular polar angle 6 normalized in the interval < 8 < n/2 [ the angles 6 > ir/2 
do not bring about anything new as f(n — 6) = f{6) ] is P{9) = svo.0. After an elementary 
transformation, we get a normalized probability in terms of /: 

P(f)df = ( 6 .4) 

As earlier, < / >= 1/3, but the distribution in / is now wide and the values / = and / = 1 
are quite probable. 

Thus, a hope exists that in, experiments with heavy ion collisions at RHIC, wild fluctuations 
in the fractions of neutral and charged pions would be observed. That would be a direct 
experimental indication that a quasi-phase-transition occurs where domains of disoriented 
chiral condensate of noticeable size are developed in a cooling stage. One can recall in this 
respect mysterious Centauro events with anomalously large fraction of neutral or of charged 
particles observed in cosmic ray experiments [[37j. Who knows, may be that was the first 
experimental observation of the QCD phase transition. 
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Fig. 5. Domains of disoriented chiral condensate. 




FIG, 2 . Presence and absence of the finite-temperature 
QCD phase transition as a function of m Uid a and m+G. Mass 
values for which the transition is and is not seen on a 16 3 x4 
lattice arc denoted respectively by solid cErclcs and squares. 
The physical point, Indicated roughly by the dashed circle, lies 
in the region of no transition. 



